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Ph 1 abstract 

Using numerical simulations based on solving the general relativistic hydrodynamic 
. equations with the CoCoNuT code, we study the dynamics of a phase transition in 

the dense core of isolated rotating neutron stars, triggered by the back bending insta- 
,— bility reached via angular momentum loss. In particular, we investigate the dynamics 
p<^ of a migration from an unstable configuration into a stable one, which leads to a 

' mini-collapse of the neutron star and excites sizeable pulsations in its bulk until it 

acquires a new stable equilibrium state. We consider two equations of state which 

■ exhibit softening at high densities, a simple analytic one with a mixed hadron-quark 
I*', phase (where the hadron pressure is approximated by a polytrope) in an intermediate 

O ■ pressure interval and pure quark matter at very high densities, and a microphysical 

l-j ' one that has a first-order phase transition at constant pressure with a jump in density, 

r/i t originating from kaon condensation. Although the marginally stable initial models are 

i i ■ rigidly rotating, we observe that during the collapse (albeit little) differential rotation 

is created. We analyze the emission of gravitational radiation in such an event, which 

, in some models is amplified by mode resonance effects, and assess its prospective de- 

■ tectability by current and future interferometric detectors. We expect that the most 
| favorable conditions for dynamic migration exist in very young magnctars. The rate 

t-H , of such events in our Galaxy is of the order of one per century and rises to about one 

per year if the Virgo cluster of galaxies is considered. We find that the damping of 

\ the post-migration pulsations and, accordingly, of the gravitational wave signal ampli- 

» . tude strongly depends on the character of the equation of state softening (either via a 

' density jump or continuous through a mixed state). The damping of pulsations in the 

Q\ . models with the microphysical equation of state is caused by dissipation associated 

<^ ■ with matter flowing through the density jump at the edge of the dense core. If at 

| work, this mechanism dominates over all other types of dissipation, like bulk viscosity 

• i-h , in the exotic-phase core, gravitational radiation damping, or numerical viscosity. 

, 

5_j ■ Key words: hydrodynamics - relativity - methods: numerical - stars: neutron - 

C$ ' stars: pulsations - stars: phase-transition - stars: rotation 



1 INTRODUCTION 

One of the mysteries of neutron stars is the structure of 
their cores where the density exceeds normal nuclear den- 
sity p nuc ~ 3 x 10 14 g cm -3 . For p < p nuc the constituents of 
dense matter are well known (neutrons, protons, and elec- 
trons). Here the respective physical theories can be con- 
fronted with nuclear physics data and are therefore under 
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control (see e.g. iHaenseh Potekhin fc YakovlevlkoOit ). How- 
ever, for p significantly larger than p nU c, the predictions of 
different plausible theories diverge with increasing density. 
The structure of matter at 5 - 10 p nU c, i.e. at densities charac- 
teristic for the center of massive neutron stars, is unknown. 
Still, in order to construct neutron star models, the equa- 
tion of state (EoS) of matter in this regime is needed as an 
input. Some theories of dense matter predict a strong soft- 
ening of the EoS at p > 2 p nuc , caused by the appearance 
of various exotic phases of matter, such as pion condensate, 
kao n condensate, or deconfined quark matter (for a review 
seelHaensel. Potekhin fc Yakovlevll2007l ; lGlendennin"3l2000l : 
IWeberlll999h . The prediction of observable signatures of the 
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presence of exotic phases in the dense core of neutron stars is 
thus crucial for testing physical theories of hadronic matter 
at extreme densities. 

A softening of the EoS is expected to influence the ro- 
tational properties of neutron stars. Today, we know some 
1800 radio pulsars (i.e. rotating magnetized neutron stars), 
which are spinning down due to angular momentum loss via 
radiation. They are observed using radio telescopes, where 
the evolution of the rotation period with time is measured, 
and some of them are also being monitored by the newly 
constructed gravitational wave detectors (albeit without a 
reliable detection so far in that case). If the EoS of matter 
in a pulsar core actually softens at some density threshold, 
this could have a dramatic impact on the structure of the 
neutron star that completely alters the pulsar timing prop- 
erties: One of the possibilities is that the pulsar could enter 
a back bending episode, which is reflected by the reversal 
of the spin-down into a s pin-up (|Glendenning. Pei fc Weberl 
ll997l : IZdunik at al.ll2006l ). If the softening of the EoS is suffi- 
ciently strong, the back bending could lead to an instability 
at which the rotating neutron s tar dynamically mi grates to 
a new stable equilibrium state ijZdunik at al.ll2006[ ). 

Depending on the mass of the neutron star, this migra- 
tion process can be accompanied by a considerable increase 
in central compactness and spin-up. Immediately after mi- 
gration, the spinning neutron star's bulk is expected to pul- 
sate with sizeable amplitude, which would result in a burst of 
gravitational radiation. This prospect gives additional mo- 
tivation for studying the dynamical migration phenomenon. 
A scenario with similar dynamics (i.e. a dynamic migration 
of rotating neutron stars towards a new stable equilibrium 
state at higher central density along with a spin-up, the exci- 
tation of strong pulsations, and the e mission of g r avitat ional 
waves) was recently inve stigated by iLin et alT (|2006l ) and 
lAbdikamalov et all (|2009h . However, while in that work the 
original, unperturbed initial configuration is in stable equi- 
librium and the transient is initated by an instantaneous 
change in the EoS (resulting in a somewhat artificial deple- 
tion of pressure at high densities) , the migration process due 
to back bending need not be initiated by an ad-hoc modi- 
fication of the EoS, as this instability starts at an unstable 
equilibrium state. 

In this work we investigate the dynamics of such a mi- 
gration of pulsar models, triggered by the instability de- 
scribed above. This is accomplished by performing fully 
nonlinear time-dependent numerical simulations with a two- 
dimensional version of the general relativistic hydrodynam- 
ics code CoCoNuT. From a multitude of existing EoSs we 
select two types which feature softening at high densities: 
The fi rst is the simple analytic MUn EoS of IZdunik at al.l 
(|2006l ) which obeys a polytropic relation at low densities, 
describes a mixed hadron-quark phase at intermediate den- 
sities (where the hadron pressure is again approximated by a 
polytrope, however with a different adiabatic exponent), and 
has pure quark matter at very high densities. The second one 
is a tabulated microphysical EoS that exhibits a first-order 
phase transition at constant pressure with a jump in density, 
originating from kaon condensation. For these two EoSs, as 
initial models we select marginally stable rotating configura- 
tions with different baryon masses, and follow the migration 
and the subsequent ring-down phase for many dynamic time 
scales. We analyze the post-migration pulsations in terms 



of amplitude, spectrum and possible damping mechanisms. 
Furthermore, we calculate the associated gravitational wave 
emission and discuss the prospects for a possible detection 
by current and future interferometric detectors like LIGO, 
VIRGO or Advanced LIGO. 

This article is organized as follows: In Section [2] we 
present the mathematical framework and numerical setup 
for the simulations, while in Section|3]we give details about 
the two EoSs employed and introduce the marginally stable 
models that undergo migration. Section U contains a survey 
of the results, with a general description of the migration dy- 
namics, an analysis of the pulsations excited by this process, 
information about the final equilibrium states, followed by 
a discussion of our models' gravitational wave emission and 
the detectability prospects. In Section [S] we summarize our 
findings and close with some conclusions. We also add some 
information about non-equilibrium effects in the EoSs in Ap- 
pendix|X]and explain the damping mechanism for pulsations 
in models with the microphysical EoS in Appendix iBl 

Unless otherwise noted, we choose dimensionless units 
for all physical quantities by setting the speed of light and 
the gravitational constant to one, c = G = 1. Latin indices 
run from 1 to 3, Greek indices from 1 to 4. 



2 FORMULATION AND NUMERICAL 
METHODS 

We construct the initial equilibrium models of rotat- 
ing neutron stars using a variant of the self-consistent 
field m ethod described in iKomatsu. Eriguchi fc Hachisul 
(1989 a|) (KEH hereafter) , as im plemented in the code RNS 
( Stergioulas fc Friedman! Il995l ) . This code solves the gen- 
eral relativistic hydrostationary equations for rotating mat- 
ter distributions whose pressure obeys a barotropic EoS. The 
resulting equilibrium models are taken as initial data for the 
evolution code. 

The time dependent numerical simulations are 
performed with the code CoCoN uT developed by 
iDimmelmeier, Font fc Muller! (|2002al lbh with a met- 
ric solver based on spe ctral methods as described in 
IDimmelmeier et al.l (|2005h . The code solves the general 
relativistic field equations for a curved spacetime in the 

3 + 1 split under the assumption of the conformal flatness 
condition (CFC) for the three-metric. The hydrodynamics 
equations are formulated in conservation form, and are 
evolved by high-resolution shock-capturing schemes. 

In the following, we present the mathematical formula- 
tion of the metric and hydrodynamics equations, and then 
summarize the numerical methods used for solving them. 
Note that this part is pr actically identical to Section 2 in 
lAbdikamalov et al.l (|2009h . 



2.1 Metric equations 

We adopt the AD M 3 + 1 formalism by 

lArnowitt. Deser fc Misnerl (|l962l ') to foliate a spacetime 
endowed with a metric g^ v into a set of non-intersecting 
spacelike hyp ersurf aces. The line element then reads 

ds 2 = g^ v dx^ dx v = —a 2 dt 2 +-y i j(dx l +(3 1 dt)(dx J +/3 J dt), (1) 
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where a is the lapse function, /3 1 is the spacelike shift three- 
vector, and 7^ is the spatial three-metric. 

In the 3 + 1 formalism, the Einstein equations are split 
into evolution equations for the three-metric jij and the ex- 
trinsic curvature Kij, and constraint equations (the Hamil- 
tonian and momentum constraints) which must be fulfilled 
at every spacelike hypersurface. 

The fluid is specified by the rest-mass density p, the 
four- velocity u M , and the pressure P, with the specific en- 
thalpy defined as h — 1 + e + P/p, where e is the specific 
internal energy. The three-velocity of the fluid as measured 
by an Eulerian observer is given by u l = u'/(au°) + /3*/a, 
and the Lorentz factor W — au° satisfies the relation 

W = l/Vl-ViV\ 

Based on the ideas of llsenberd (|2008l ) 



and 



IWilson. Mathews fe Marronettil (Il996), and a s it wa s done 
in the work of iDimmelmeier. Font fc Mullerl (|2002al lbh. we 



approximate the general metric g^ v by replacing the spatial 
three-metric 7^ with the conformally flat three-metric 



7<i = 9 Jij, 



(2) 



where jij is the flat metric and <j> is a conformal factor. In 
this CFC approximation, the ADM equations for the space- 
time metric reduce to a set of five coupled elliptic non-linear 
equations for the metric components, 



A(j> = -2tt0 5 (phW 2 - P) 



7K ■ 

A(a<t>) = 2wact> 5 (ph(3W 2 - 2) + 5P) + acj> 5 ' J , (3) 



A/3* = l6-Ka<t>phW 2 v + I^K^Vj (aq 



1 



V 8 V fc /? fc 



where the maximal slicing condition, K\ = 0, is imposed. 
Here V< and A are the flat space Nabla and Laplace oper- 
ators, respectively. For the extrinsic curvature we have the 
expression 



K, 



277 



(4) 



which closes the system ©. 

We rewrite the above metric equations in a mathemati- 
cally equivalent form by int roducing an auxiliary vector field 
W i and obtain ||Saiicll2004h 



A</> = -2tt0 5 (phW 2 - P) 



A(a(j>) = 2na(j> 5 (ph(3W 2 - 2) + 5P) + a<j> 



■7K.K- 
8 ' 



A/3 4 = 2V : 



(5) 



6 K lJ ) - ^V l V fe /f 



(6) 



AW 1 = 8TT(f> w phW 2 v l - -X7 l X7 k W k , 
where the flat space extrinsic curvature is given by 

kij = ViWj + %Wi - ^fijV k w k 

and relates to the regular extrinsic curvature as Kij = 
4> 2 K i:j and K ij = (t> 10 K i3 . The virtue of this reformu- 
lation of the metric equat ions is discussed in detail by 
I Codero- Carrion et all (|2009h . 



Note that the metric equations do not contain explicit 
time derivatives, and thus the metric is calculated by a fully 
constrained approach, at the cost of neglecting some evo- 
lutionary degrees of freedom in the spacetime metric (e.g. 
dynamic gravitational wave degrees of freedom). 

The accuracy of the CFC approximation has been 
tested in various works, both in the context of stellar 
core collapse and for equilibrium models of neutron stars 
(for a detailed comparison betw een the CFC approxima- 
tion and full gen eral relativity, see IDimmelmeier et al.ll2006l ; 
lOtt et aT1l2007l . and references therein). The spacetime of 
rapidly (uniformly or differentially) rotating neutron star 
models is still very well approximated by the CFC met- 
ric (J2Jl. The accuracy of the approximation is expected to 
degrade only in extreme cases, such as a rapidly rotating 
black hole or compact binary systems. 

2.2 General relativistic hydrodynamics 

The hydrodynamic evolution of a relativistic perfect fluid 
is determined by a system of local conservation equations, 
which read 



0. 



0. 



where J M = pit M is the rest-mass current, and V,j denotes 
the covari ant derivative with re spect to the four-metric <; M „. 
Following iBanvuls et Tall (|l997l ) we introduce a set of con- 
served variables in terms of the primitive (physical) variables 
(fi,Vi,e): 

D = pW, S l =phW 2 v 1 , T = phW 2 -P-D. (8) 

Using the above variables, the local conservation laws Q 
can be written as a first-order, flux-conservative hyperbolic 
system of equations, 



-gs, 



dt dx 1 
with the state vector, flux vector, and source vector 

U = [D,S j ,r], 

F l = \pv l , SjV + 5}P, tv* + Pv 



(9) 



'2 dxi 



+ 



(10) 



T 



x \fli with 



respectively. Here v 1 — v l — (3 1 /a, and \/~^g 
g = det(g My ) and 7 = det(7y). 

The system of hydrodynamics equations <(9j is closed by 
an EoS, which relates the pressure to some thermodynami- 
cally independent quantities, in our case P = P(p,e). 

2.3 Numerical methods for solving the metric and 
hydrodynamics equations 

The hydrodynamic solver performs the numerical time inte- 
gration of the system of conservation equations (|9} using a 
high-resolution shock-capturing (HRSC) scheme on a finite- 
difference grid. In (upwind) HRSC methods a Riemann 
problem has to be solved at each cell interface, which re- 
quires the reconstruction of the primitive variables (p, v l , e) 
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at these interfaces. We use the piecewise parabolic method 
(PPM) method for the reconstruction, which yields third- 
order accuracy in space for smooth flows and away from 
extrema. The numerical fluxes are co mputed by means o f 
Marquina's approximate flux formula (|Donat et al.|[l998h . 
The time update of the conserved vector U is done using 
the method of lines in combination with a Runge-Kutta 
scheme with second-order accuracy in time. Once the state 
vector is updated in time, the primitive variables are re- 
cover through an iterative Newton-Raphson method. To 
numerically solve the elliptic CFC metric equations (0 we 
make use of an iterative non-linear solver based on spectral 
methods. The combination of HRSC methods for the hy- 
drodynamics and spectral methods for the metric equations 
in a multidimensional num erical code has been in detail in 
iDimmelmeier et al.l (120051 ). 

The CoCoNuT code utilizes Eulerian spherical polar 
coordinates {r, 8}, as for the models discussed in this work 
we assume axisymmetry with respect to the rotation axis 
and additionally equatorial symmetry. The finite-difference 
grid consists of 160 radial and 40 angular grid points, which 
are equidistantly spaced. A small part of the grid, which 
initially corresponds to 20 radial grid points, covers an ar- 
tificial low-density atmosphere extending beyond the stellar 
surface, whose rest-mass density is 10~ 17 of the initial cen- 
tral rest-mass density of the star. 

Since the calculation of the spacetime metric is 
computationally expensive, the metric is updated only 
once every 25 hydrodynamic time steps during the evo- 
lution and extrapolated when needed. The suitability 
of this procedure is tested and d iscussed in detail in 
IDimmelmeier. Font fc Miillerl (|2002al ). We also note that 
tests with different grid resolutions were performed to as- 
certain that the regular grid resolution specified above is 
appropriate for our simulations. 



2.4 Gravitational waves 

The gravitational waves emitted by the collapsing neutron 
star are computed using the Newtonian quadrupole for- 
mula in its first time-integrated form (the first-moment 
of momentum density formulat i on as described in detail 
in IDimmelmeier. Font fc Miillerl l2002bh in the variant of 
IShibata fc Sekiguchil (|2004f ). It yields the quadrupole wave 
amplitude A^q as the lowest order term in a multipole ex- 
pan sion of the rad iation field into pure-spin tensor harmon- 
ics (|Thornelll980l ). The wave amplitude is related to the 
dimen sionless gravitational wave strain h i n the equatorial 
plane (|Dimmelmeier. Font fc Miillerl l2002bh 



1 1 15 A™ 



h = —\ — 

8V 7T 



8.8524 x 10" 



Af 2 lOkpc 
10 3 cm r 



(11) 



with r being the distance to the emitting source. 

We point out that although the quadrupole formula is 
not gauge invariant and is only strictly valid in the Newto- 
nian slow-motion limit, for gravitational waves emitted by 
pulsations of rotating neutron stars it yields results that 
agree very well in phase and to about 10% -20% in ampli- 
tude with more sophistic ated methods (|Shibata fc Sekiguchil 
120031 ; iNaear et alJl2007T ) . 



3 MODEL SETUP 



3.1 Equations of state 



One of our two sets of models is based on a simple analytic 
three-phase EoS which is identical to the m ixed unstable 
(MUn) EoS discussed in lZdunik at~aH (|2006l "FI: 



In the normal hadronic phase, i.e. for p < pi = 3.32 x 
(corresponding to n\ = 0.2 x 10 39 cm~ 3 ), the 
EoS is given by a polytrope, 



10 14 g cm -3 



Kp 1 



(12) 



with 7 = 7h = 2.5 and a polytropic constant K = — 
1.0506 x 10~ 2 (in cgs units, which equals 0.025 in units 
of pc 2 /n 7h wi th p = 1.66 x 10 14 g cm" 3 and n = 0.1 x 
10 39 cm" 3 , see IZdunik at all I2006T ). The internal energy is 
calculated according to the ideal gas EoS as 

(13) 



P(7-1) 

• In the mixed phase, i.e. for pi ^ p < P2 
10 14 g cm" 3 (corresponding to n 2 = 0.65 x 10 39 cm 



10.79 x 
the 



-3\ 



EoS is again a polytrope, now with 7 = 7 m = 1.3 and 
a po lytropic constant K = K m = 2.7976 x 10 15 (|Zdunikl 
l2000h . which ensures continuity of P at pi. The expression 
for the internal energy is more complicated than in the reg- 
ular phase, 

£=—, tt H 1, (14) 

p(7-l) m u 

where the mean atomic mass 

Pi 7h - 7m 



Pi (7h - l)(7m - 1) 



(15) 



in the mixed phase is determined by demanding that also e 
is continuous at pi. Here Pi is the pressure at pi. 

• At high densities p p2 , we assume a pure quark mat- 
ter phase, wh ose pressure is given by the MIT bag model 
|Zdunikll20o"oT ), 



P= -(£-£ 2 ) + P 2 , 



P2 



4P + £ 2 - 3P 2 
£2 + Pi 



3/4 



(16) 



where £ — p(l + e) is the total energy density of the quark 
matteiQ and P2 and £ 2 are the pressure and energy density 
at p 2 , respectively, which again enforces continuity of the 
EoS at the transition point. 

The other set of models uses an microphysical EoS, 
which exhibits a phase transition to a kaon condensate ac- 
companied by a density discontinuity. The composition of 
this EoS is as follows: 

• For the description of the inner crust (nuclei immersed 
in neutron and electron gas) and the outer crust (nuclei im- 
mersed in electron gas) at densities 5 x 10 10 g cm" 3 ^ p ^ 



1 Note that here we give the EoS in terms of the rest-mass den- 
sity p (the quantity used by our evolution code) rather than the 
baryon number density n = p/m n , where m u = 1.66 X 10 — 24 g is 
the baryon mass. 

2 By assuming that e = £ / p — 1 we can continue to use p and e 
in our system of hydrodyna mic equations for a m ore detailed 
discussion of this issue, see lAbdikamal ov et al. I <2009l) . 
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Figure 1. Dependence of the pressure P on the density p for the 
simple MUn EoS (solid line) and the microphysical EoS (dashed 
line). The transition densities p\ and p2 for both the MUn EoS 
and the microphysical EoS are indicated by the dotted lines. Note 
the domain of zero pressure slope in the microphysical EoS. 



1.5 x 10 14 g cm 3 we use the model of iDouchin fc Haensell 
(2001). In this model, neutron drip occurs at a den- 
sity of 3 x 10 11 g cm -3 . The outer crust at densities 
10 s g cm" 3 < p 5x 10 10 g cm" 3 is modeled by the EoS of 
lHaensel fc Pichonl fl994!). Finally, the low-dens ity envelope 
with p < 10 8 g cm" 3 is described by the EoS of lBavm et alj 
d!97lh . 

• The constituents of the moderately dense part of the 
core below the lower end of the density jump at pi = 
6.31 x 10 14 g cm" 3 are neutrons, protons, electrons, and 
muons. The nucleon component is described by a relativis- 
tic mean-field model with scalar s elf-coupling, as constructed 
bv lZimanvi fc Moszkowskil (|l990f ). The values of the meson- 
nucleon coupling constants are g a lm a = 3.122 fm, g^/m^ = 
2.1954 fm, and g p /m p — 2.1888 fm. The dimensionless co- 
efficients in the cubic and quartic terms in the scalar self- 
coupling are b = -6.418 x 10~ 3 and c = 2.968 x 10~ 3 , re- 
spectively. 

• Although ideally the pressure P should remain constant 
across the density jump from pi to p2 = 10.24 x 10 14 g cm" 3 , 
in order to avoid numerical problems with the hydrodynam- 
ics scheme of the CoCoNuT code we introduce a small linear 
pressure increase of 10% in this density interval. 

• For densities higher than the upper end of the density 
jump at p2 we consider kaon-condensed matter. The cou- 
pli ng of kaons to nucleons is done a c cordin g to the model 
of iGlendenning fc Schaffner-Bielichl (|l999l ) with U^ n = 
-110 MeV. 

The properties of these two EoSs are exemplified in 
Fig. [T] where we plot the dependence of the pressure P on 
the rest-mass density p. In the case of the microphysical EoS 
the approximately constant pressure between pi and p2 is 
apparent. In the dense phase, both EoSs have a very similar 
shape. 

Our approach assumes that the fluid is in thermody- 
namic equilibrium everywhere, even during the dynamic mi- 



gration process. The justification of this assumption is given 
in Appendix IA1 

As both the simple MUn EoS and the microphyical EoS 
are barotropic, there is no need to evolve the quantity r in 
the hyperbolic set of hydrodynamic evolution equations (j9] 
UOp . During the recovery of the primitive variables, the as- 
sociated quantity e is simply evaluated from the EoS using 
the rest-mass density p, and then r is computed straightfor- 
wardly according to Eq. JH}. 



3.2 Pre-migration initial models 

As shown by IZdunik at all (|2006h . for the MUn EoS se- 
quences of uniformly rotating neutron star models with con- 
stant baryon mass Mo exhibit a back bending phenomenon. 
As shown in Fig. [5] when plotting the total angular momen- 
tum J against the rotation frequency /i, for each sequence 
there exists an unstable segment between a local minimum 
and maximum of J (marked by dashed lines). As an iso- 
lated neutron star spins down at constant Mo along this 
evolutionary path, it will migrate from the last marginally 
stable configuration (denoted by a circle) to the location on 
the stable branch with the same total angular momentum 
J, provided that the migration proceeds sufficiently rapid 
that J is conserved and that the migration does not in- 
troduce a deviation of the rotation profile from uniformity. 
While the first assumption is justified, as the migration is 
a fast dynamic process, the condition of uniform rotation 
after the migration has to be verified by numerical simula- 
tions. If the rotation profile stays approximately uniform, 
then the migration path will proceed closely along the dot- 
ted lines in Fig. [2] For the microphysical EoS, we find the 
same back bending phenomenon for sequences of constant 
Mq, as presented for a slightly different selection of values 
for the gravitational mass in Fig. [3] 

For our simulation of the dynamic migration due to back 
bending, we construct initial models that correspond to the 
last marginally stable configuration marked by circle^] in 
Figs. [2] and [3] The properties of these initial models, which 
utilize the simple MUn EoS (family US) or the microphysical 
EoS (family UM) are summarized in Table [T] Following the 
assumption that during the spin down of isolated neutron 
stars viscous processes drive the rotation profile towards uni- 
formity, we initially assume rigid rotation. Then each model 
is uniquely specified by its central energy (or alternatively 
rest- mass) density and its rotation rate Ti/|Wi|, which is the 
ratio of rotational energy to gravitational binding energy. 

A remarkable property of the marginally stable initial 
models with the microphysical EoS is that their central den- 
sity is constant irrespective of the baryon mass (see Table[T|. 
This is a consequence of the first-order phase transition in 
this EoS with a density jump between pi and p2- This tran- 
sition destabilizes the neutron star exactly when a (theoret- 
ically infinitesimally small) core of the dense phase begins 



3 Even though Fig. suggests that for models with the micro- 
physical EoS, the local minimum of J for the marginally stable 
models coincides with the local minimum of /;, a magnification 
of this plot reveals that this is not the case and the J—f\ curves 
actually are smooth everywhere. 
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Figure 2. Dependence of the total angular momentum J on the 
rotation frequency /; for uniformly rotating equilibrium models 
with the simple MUn EoS for various sequences of constant total 
baryon mass Mo. Solid/dashed lines specify the stable/unstable 
segments of the curve. The locations of the marginally stable ini- 
tial models of the family US are marked by circles, and their tran- 
sition path to the stable segment assuming conservation of J and 
uniform rotation, resulting in a change A/ in rotation frequency, 
is indicated by dotted lines. 

to form in its centeiQ. In contrast, the MUn EoS does not 
posses a first-order phase transition but rather a continuous 
transition through a mixed phase without density jump. Al- 
though the mixed phase also softens the EoS, this softening 
is not as abrupt as in the microphysical EoS, but "propor- 
tional" to the size of the soft core of the mixed phase. Thus 
the star loses stability only if such a core of sizeable mass 
is already present in the center, and the exact value of p c 
when the marginally stable configuration is reached depends 
(albeit only weakly) on the mass of the neutron star. 

If the hydrodynamic and metric equations could be 
solved to arbitrarily high precision, then the marginally sta- 
ble models would already migrate if subjected to an in- 
finitely small perturbation. In our numerical setup, however, 
truncation and round-off errors together with the small er- 
rors from using the CFC approximation of the exact metric 
equations introduce effects which make it hard to predict 
if models in the close vicinity of a stable configuration are 
driven over the edge of stability or not. In order to ensure 
that in our numerical simulations the initial models undergo 
the migration robustly in a finite evolution time, we choose 
to initially apply a finite perturbation to the rest-mass den- 
sity of the form 



Ap = a sin [ — — 'i 



(17) 
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Figure 3. Same as Fig. [2] but for models with the microphys- 
ical EoS. The locations of the marginally stable initial models 
of family UM are marked by circles. At the models marked by 
plus symbols, the Kepler limit interrupts the two sequences with 
M = 1.55 Mq and M = 1.60 Mq. The transition path of the 
Mo = 1.55 Mq sequence is not affected by this discontinuance, 
while for the Mo = 1.60 Mq sequence there exists no model on 
the stable branch with the same J as the marginally stable one. 
Note that a fiducial line connecting all marginally stable models 
with this EoS would end with an Mo = 1.639 Mq model rotating 
at the Kepler limit (marked by an x symbol). 
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Figure 4. Location of the three evolved stable models (marked 
by crosses) in the J—f\ diagram for uniformly rotating equilibrium 
models with the simple MUn EoS and Mo = 1.55 Mq (see also 
Fig. EJl. 



where a = 10 3 is the perturbation amplitude and r atm is 



4 For numerical reasons, in our study the density discontinuity 
in the microphysical EoS is slightly softened by introducing a 
small pressure gradient, which leads to a threshold value p c j = 
9.8 X 10 14 g cm -3 that is a bit lower than p2, where the dense 
phase begins. The smooth dependence of J on /j discussed in 
Footnote [3] is also a consequence of this pressure gradient. 



the (latitude dependent) boundary of the neutron star to 
the surrounding low-density atmosphere. In order to assess 
the evolutionary behavior of models near the stability limit, 
we augment model US3 by the three models of family SS3 
which are located close to it still on the stable branch, as 
shown in Fig. [4] 
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Table 1. Summary of the set of migrating marginally stable initial models with the simple MUn EoS (family US) and the microphysical 
EoS (sequence UM) as well as the stable initial test models with the simple MUn EoS (family SS3). Mo is the total baryon mass, M is 
the total gravitational mass, € c ,i is the initial central energy density, p Cj i is the initial central rest-mass density, /j is the initial rotation 
frequency, J is the total angular momentum, r p ;/r c ; is the initial ratio of polar radius to equatorial radius, and / 1 Wi is the initial 
rotation rate. 



Model 


M 
[M e ] 


M 

[M ] 


S c ,i 

[10 14 g cm" 3 ] 


Pc,i 

[10 14 g cm" 3 ] 


Pc,i/Pl 


Pc,i/P2 


fi 
[kHz] 


J 

[GM^c- 1 ] 


r P,i/ r o,i 


Ti/\Wi 
[%] 


US1 


1.46 


1.35 


10.22 


9.05 


2.73 


0.84 


0.057 


0.07 


0.998 


0.04 


US2 


1.50 


1.39 


10.39 


9.19 


2.77 


0.85 


0.364 


0.51 


0.939 


1.59 


US3 


1.55 


1.44 


10.39 


9.19 


2.77 


0.85 


0.505 


0.77 


0.880 


3.19 


US4 


1.60 


1.48 


10.22 


9.05 


2.73 


0.84 


0.595 


1.00 


0.829 


4.64 


US5 


1.65 


1.53 


10.27 


9.09 


2.74 


0.84 


0.665 


1.22 


0.781 


6.03 


UM1 


1.35 


1.26 


11.03 


9.80 


2.95 


0.91 


0.340 


0.36 


0.953 


1.19 


UM2 


1.40 


1.30 


11.03 


9.80 


2.95 


0.91 


0.528 


0.61 


0.885 


2.98 


UM3 


1.45 


1.35 


11.03 


9.80 


2.95 


0.91 


0.646 


0.82 


0.823 


4.63 


UM4 


1.50 


1.39 


11.03 


9.80 


2.95 


0.91 


0.732 


1.01 


0.766 


6.16 


UM5 


1.55 


1.44 


11.03 


9.80 


2.95 


0.91 


0.798 


1.21 


0.711 


7.58 


UM6 


1.60 


1.48 


11.03 


9.80 


2.95 


0.91 


0.848 


1.39 


0.656 


8.84 


SS3a 


1.55 


1.44 


8.21 


7.39 


2.23 


0.68 


0.492 


0.82 


0.873 


3.46 


SS3b 


1.55 


1.44 


8.87 


7.94 


2.39 


0.74 


0.490 


0.79 


0.879 


3.27 


SS3c 


1.55 


1.44 


9.60 


8.54 


2.57 


0.79 


0.494 


0.77 


0.881 


3.18 



4 RESULTS 

4.1 Dynamics of the migration process 

It is well known that the migration of a marginally stable 
nonrotating spherical neutron star model with a polytropic 
EoS and adiabatic coefficient 7 = 2 leads to an expansion 
of t he star (for a nu merical simulation of this scenario, see 
e.g. lFont etai1l2002l ). On the unstable branch, such a model 
possesses a very compact structure with high central den- 
sity p c , while the corresponding model with same baryon 
mass (and same gravitational mass) on the stable branch 
has a much smaller value of p c . Consequently, the entire 
star significantly expands in size and, after some ring-down 
pulsations following the initial expansion, settles down to a 
less dense new equilibrium configuration. 

In stark contrast, in our case where two different EoSs 
are used and the marginally stable initial neutron star mod- 
els rotate rapidly, a collapse instead of an expansion is 
expected. This follows directly from the fact that during 
the migration the rotation frequency / increases as seen in 
Figs. [2] and [3] while the total angular momentum is con- 
served, if the uniform rotation profile is approximately pre- 
served throughout the evolution (which is a good approxi- 
mation, as later shown in Section [4.2|) . Still, as in the case 
of a nonrotating simple polytrope the new equilibrium state 
- now characterized by an increase in central density com- 
pared to the initial value - is overshot after the primary 
contraction and is reached only after some ring-down pulsa- 
tions. 

Before we present the collapse dynamics of the 
marginally stable models that undergo migration, we first 
discuss the evolution of models from the stable branch that 
are close to but still sufficiently far from the stability limit, 
like models SS3a or SS3b. As expected such models remain 
stable throughout the evolution and simply oscillate around 
their equilibrium state due to the small perturbation ap- 
plied initially (see Fig. [SJ. However, the amplitude of these 
damped pulsations already grows from model SS3a, which 
is the one farthest from the marginally stable model US3, 




10 20 30 

t [ms] 



Figure 5. Time evolution of the central density p c for the 
marginally stable migration model US3 (solid line) and the three 
stable models SS3a, SS3b, and SS3c (dashed, dotted, and dash- 
dotted lines, respectively) with the simple MUn EoS and Mo = 
1.55 Mq. Note that model SS3c actually migrates as well, result- 
ing from the finite initial perturbation along with errors intro- 
duced by the numerical treatment and the CFC approximation. 



to model SS3b that is closer to US3. Interestingly, model 
SS3c, which lies very near model US3 but is still stable, al- 
ready undergoes the collapse to higher central densities and 
more rapid rotation that is associated with the migration 
to the opposite stable segment. In this case the errors from 
the numerical treatment and the CFC approximation, and 
in particular the small but finite initial perturbation suffice 
to push the neutron star model over the stability limit. Its 
subsequent central density evolution is very similar to that 
of the migrating marginally stable model US3, which is also 
shown in Fig. 

When comparing the time evolution of the central den- 
sity for different models of the same family, as presented in 
Figs. [6] and [7] for the two EoSs, one recognizes a generic be- 
havior. Within slightly more than 1 ms, p c increases from 
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its initial value to a peak of p c . p ~ 15-16 x 10 g cm , 
and subsequently undergoes damped harmonic oscillations 
around an asymptotic value that lies between p c , P ~ 13 and 
14 x 10 14 g cm -3 . The actual numbers for the density in- 
crease in the peak and the final value are given in Tabled 

The baryon mass fractions in Table [2] also demonstrate 
that initially no fluid element in the neutron star has a den- 
sity above p2, which is consistent with the values for p c> ; 
given in Table [T] while around 35% and 60% of the stellar 
mass is at a density above pi for the MUn EoS and the 
microphysical EoS, respectively. Later in the evolution, the 
mass of the star above pi grows noticeably, and even a sig- 
nificant fraction of Mo has a density greater than p2. 

When comparing Figs. [6] and [7] it is apparent that while 
the density oscillations of models with the MUn EoS induced 
by the migration are somewhat damped, the damping in the 
models with the microphysical EoS is much strongeiQ After 
only a few cycles, the pulsation have died out and p c (as well 
as other hydrodynamic and metric quantities at various loca- 
tions in the star) remain practically constant. By performing 
simulations with significantly higher grid resolution, we have 
ensured that this effect is not caused by insufficient numer- 
ical accuracy. In addition, since our code utilizes Riemann 
flux solvers, which can handle discontinuities in the hydro- 
dynamic quantities without p roblem (see e.g. the relativistic 
shock tube test presented in iDimmelmeier. Font fc Mullerl 
l2002al ). we can exclude that the pressure jump in the micro- 
physical EoS leads to numerically questionable results. By 
using a relative convergence criterion of 10 -12 , we also as- 
certain that the recovery routine for the primitive quantities 
is performed to sufficient precision. 

When investigating the strong damping of post- 
migration pulsations in the models of family UM with the 
microphysical EoS, we find strong evidence that it is a re- 
sult of a very efficient physical damping mechanism that is 
reflected in our code. As shown in detail in Appendix[B] the 
kinetic energy stored in the pulsations is dissipated during 
the compression phase of each pulsation, as matter passes 
through the density jump from pi to pi associated to the 
first-order phase transition in this EoS. While in reality 
matter flows into the dense core at supersonic speed and is 
heated up by shock dissipation, in our simulations the flow 
remains subsonic. Still, the large density jump results in a 
rather steep radial inflow velocity profile, leading to numer- 
ical dissipation of comparable strength. If the density jump 
in the EoS is gradually removed, dissipation continuously 
decreases to the normal level expected for a code based on 
HRSC methods. 

Furthermore, non-equilibrium processes produce en- 
tropy via heating. This can be represented by bulk viscos- 
ity, but as the corresponding additional dissipation is much 
smaller than the dissipation caused by the phase transition 
(see Appendix IB3|I . it is admissible to neglect this mecha- 
nism as it is done in our simulations. Another, less impor- 
tant mechanism that is responsible for the damping of post- 
migration pulsations is the continuous but weak conversion 
of kinetic energy stored in the pulsations into internal energy 
by numerical dissipation in the bulk of the star (excluding 
the parts of the star that are affected by a possible density 



5 Note the different time scales in the two figures! 



jump in the EoS). However, as this process is equally active 
in all models irrespective of the EoS, in the models with the 
microphysical EoS it can at most explain damping in a small 
extent similar to what is observed in the models with the 
simple MUn EoS, which have a much larger damping time 
scale. In contrast to som e of the rapidly rotating neutron 
star models discussed by lAbdikamalov et al.l (|2009h . which 
undergo a phase-transition-induced collapse and whose post- 
collapse pulsations are mainly damped by the shedding of 
mass from the surface to the atmosphere, we do not find 
any evidence for the effect of damping by shedding in any 
of our simulations. The loss of mass across the boundary 
of our neutron star models remains well below 1% until we 
deliberately terminate the evolution at ti — 50 ms. 

It is remarkable that according to Table [2] all respective 
models of family US and UM show a generic "radial" col- 
lapse behavior, expressed by a very similar relative increase 
in central density and fractions of baryon mass above the 
respective (EoS dependent) transition densities at the col- 
lapse peak and the final equilibrium, although the change 
A/ in rotation frequency along the migration path in the 
J—fi diagram is rather different (see Figs. [2] and [3} . For the 
two investigated EoSs, the balance between pressure, gravi- 
tational, and centrifugal forces already lead to almost identi- 
cal values of these mass fractions and the central densitjjf] in 
the marginally stable initial models (see also Table [TJ, and 
apparently this does not change during the collapse. We also 
point out that while the migration lets p c increase by about 
60 - 70% to its peak value and 40 - 50% to the final equilib- 
rium value, the relative changes in the rotation rate Ti/|Wi| 
are comparatively small with values of about 9-12% and 
6-7%, respectively. 

4.2 Final equilibrium states after the migration 
and spectrum of the pulsations 

If the uniformity of the profile rotation is maintained dur- 
ing the dynamic migration and thus the collapsed neutron 
star in its new final e quilibrium state also rotates rigidly, as 
already speculated bv lZdunik at al.l (|200rj ). then the migra- 
tion essentially follows the path marked by the horizontal 
dotted lines in Figs. [5] and [3] modified only negligibly by a 
tiny change in J due to gravitational wave emission (which 
however cannot be modeled with our formulation of the met- 
ric equations). 

With our fully dynamic code we are now able to test this 
assumption of preservation of the rotation profile during the 
collapse and the subsequent ring-down. In Figs. [5] and [5] the 
radial profiles of the rotation velocity v v in the equatorial 
plane are plotted both for the initial model (where they 
obviously scale linearly with radius r) and at the end of 
the evolution at 50 ms for the simple MUn EoS and the 
microphysical EoS, respectively. In both cases the degree of 
differential rotation created by the migration is rather small, 
although the central density increases by up to 50%. It is 
actually possible to approximate the final rotation profile 
rather accurately by an analytic law for the specific angula r 
momentum j (see e.g. lKomatsu. Eriguchi fc Hachisull989al ). 

6 In the case of the microphysical EoS, p c i is indeed independent 
of Mq as explained in Section l3~2l 
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Table 2. Summary of the migration dynamics for models with the simple MUn EoS (family US) and the microphysical EoS (sequence 
UM). p CjP is the peak central density obtained at the first collapse peak, while p c f is the asymptotic value for the central density after 
the pulsations have died out and p CjS b is the central density of the corresponding equilibrium model with the same value of J from 
the stable branch in Fig, l2l and |31 respectively (for model UM6 that corresponding model is beyond the Kepler limit). Mq p>Pi ^ 2 ^/Mq, 
A^0,p>p 1 y 2 .p/^0> an< ^ ^0,p>p 1 / 2 .f/^° are factions of baryon mass above the (EoS dependent) transition densities pi and pi at the 
initial time, the time of the maximum central density peak, and averaged at late evolution times, respectively. |/i| max is the maximum 
gravitational wave strain, with the values without taking into account mode resonance effects given in parentheses. 
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0.35 


0.68 


0.57 
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0.33 
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0.34 
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0.54 


0.26 
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Figure 6. Time evolution of the central density p c for three migration models of family US with the simple MUn EoS. 
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Figure 7. Time evolution of the central density p c for three migration models of family UM with the microphysical EoS. 
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Figure 8. Radial profiles in the equatorial plane of the rotation velocity for three models with the simple MUn EoS. The solid and 
dashed lines are at the initial and final evolution time, respectively, for migration models of family US, while the dotted lines show 
at the initial evolution time for stable post-migration models of family PMS. 
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Figure 9. Radial profiles in the equatorial plane of the rotation velocity for three models with the microphysical EoS. The solid and 
dashed lines are at the initial and final evolution time, respectively, for migration models of family UM, while the dotted lines show 
at the initial evolution time for stable post-migration models of family PMM. 



j = A 2 (n c -n), (is) 

where A parametrizes the degree of differential rotation 
(stronger differentiality with decreasing ^4) and Q c is the 
value of the angular velocity il at the center. This simple 
relation is widely used in general relativistic hydrodynamics 
to model differential rotation in compact stars. The differen- 
tial rotation length scale A can be put into relation with the 
equatorial radius r c in isotropic coordinates by defining the 
dimensionless quantity A. If A <C 1, then the rotation profile 
is very differential (with a larger specific angular momentum 
distribution at small distances to the rotation axis), while 
i > 1 leads towards the limit of uniform rotation. For the 
approximated post-migration rotation curves of our models 
we find values A that are always larger than 1 and reach 
almost 3 for some models (as presented in Table [3} which 
confirms the conclusions from the visual analysis of Figs. [8] 
and[9] Despite the creation of some differential rotation dur- 
ing the dynamic migration, the post-migration models are 
thus rather well approximated by the equilibrium models on 
the stable branch in the J-fi diagram (given by the right end 



points of the dotted lines in Figs.[2]and[3} that have the same 
value of J as the marginally stable initial model. Further ev- 
idence for this conclusion is given by the observation that 
the final central density p Cj f of any post-migration model 
differs from the central density p c ,sb of the corresponding 
model on the stable branch by at most 3% (see Table [2}. 

Although the pulsations triggered by the initial collapse 
are only slowly damped for the models with the simple MUn 
EoS, they will eventually die out, and a stable equilibrium 
will be reached. With the knowledge of an approximate an- 
alytic solution for the final rotation profile we are thus able 
to compute stable equilibrium models which should closely 
match the post-migration configurations obtained from the 
dynamic evolution. For this we utilize the same method 
employed to construct the marginally stable initial models 
for the migration simulations, with the important difference 
that now slightly differential rotation profiles with the values 
for A from Tableware needed. Using the same EoS as for the 
respective migration model, only one additional parameter 
to uniquely determine each of these new equilibrium models 
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is left, like e.g. the baryon mass Mo, the gravitational mass 
M, the total angular momentum J, or the central (energy 
or rest-mass) density £ c or p c . 

For practical reasons we chose p c as the parameter of 
choice, which we select such that it equals the value of p c ,f 
in Table [2] With this method for each migration model of 
the two families US and UM an accompanying equilibrium 
model that resembles the final state is constructed. We label 
the resulting families of approximative post-merger models 
by PMS and PMM for the simple MUn EoS and the micro- 
physical EoS, respectively. The parameters for these mod- 
els are summarized in Table [3] The excellent agreement in 
both Mo and M as well as J between the respective mem- 
bers of families of marginally stable initial models and post- 
migration equilibrium models (compare the values given in 
Tables [1] and [3} illustrates the validity of this approach^. 
Therefore, it is possible to very accurately describe the neu- 
tron star at the end point of the migration by a simple equi- 
librium model obeying an analytic rotation law, which fa- 
cilitates a variety of ways to further analyze the outcome of 
the migration process. 

This method even works for the peculiar model UM6 
with the microphysical EoS, whose marginally stable initial 
configuration has no counterpart on the stable branch at 
constant J and respectively higher / in the J-fi diagram for 
uniformly rotating equilibrium models in Fig.|3j as models in 
this segment of the sequence are beyond the Kepler limit. It 
is noteworthy that model UM5 exhibits the same migration 
dynamics as the other models of its family, without contract- 
ing to excessively high densities (or even collapsing a black 
hole) or being subject to mass shedding. The small (albeit 
nonzero) development of a differential rotation profile dur- 
ing the migration suffices to create a stable final equilibrium 
state that is located within the Kepler limit. We are thus 
able to construct the corresponding post-migration equilib- 
rium model PMM6, whose value A = 2.5 for the differential 
rotation parameter is typical for family UM (see Table 

We find that even the extremal marginally stable model 
with Mo = 1.639 Mg, which already initially rotates at 
the Kepler limit (see Fig. [3} migrates as the other models 
to a stable, differentially rotating equilibrium configuration 
without mass shedding. The segments of the Mo = const, 
sequences in the J-fi diagram that are beyond the Kepler 
limit are therefore irrelevant for the actual migration pro- 
cess, as the fiducial migration path indicated by dotted lines 
in that figure requires exact preservation of uniform rotation 
throughout the migration, which is not obeyed during a fully 
dynamic evolution. We note that the marginally stable mod- 
els with the simple MUn EoS also end at the Kepler limit 
with a model that, however, has a rather large baryon mass 
Mo = 1.899M0. As this model's total angular momentum 
J = 2.30 is very large as well, we refrain from including it 
in the respective J-fi diagram in Fig. [2] 



4.3 Analysis of the pulsations modes 

A clear example of the benefit of being able to approxi- 
mate the final post-migration states of our models by sim- 
ple equilibrium configurations is the possibility to compre- 
hensively examine the quasi-normal pulsation modes trig- 
gered by the migration by perturbing those respective post- 
migration equilibrium models. It is difficult in some cases to 
unambiguously identify the mode nature and its eigenfunc- 
tion structure directly from the dynamic simulation of the 
migration itself, as both the absolute and relative excitation 
of the various modes is not under control, and the pulsation 
spectrum is rather complex (see the solid lines in Fig. llOp . In 
contrast, in the corresponding equilibrium models, by apply- 
ing a initially a selection of small trial eigenfunctions of de- 
fined parity, we can excite a range of particular quasi-normal 
modes in the linear regime with controlled amplitude. For 
I = oscillations, we set the radial velocity to 



v r = a sm | 7r 



(19) 



with amplitude a = —0.005, and for I = 2 we set the angular 
velocity to 



(20) 



vg = a smtf cosy sm | -k 



with a = 0.01. In order to improve the results obtained 
from theseperturbed equilibrium models even more, in dif- 
ficult caseqj we apply a further recycling run for each of the 
main quasi-normal modes, which are the fundamental quasi- 
radial F-mode, the fundamental quadrupolar 2 /-mode, and 
their respective harmonics, the .Hi-mode and the 2 pi-mode 
(for this we essentially follow the tec hnique described in 
iDimmelmeier. Stergioulas fc F ont 2006). This method facil- 
itates an accurate extraction of both the frequency and the 
eigenfunction of each mode, and thus allows us to clearly 
identify the mode nature. The improvement in the clarity of 
the pulsation spectrum from using the trial eigenfunctions 
alone is apparent from the dashed lines in Fig. 1101 where 
1 — 0/2 parity was used in the left/right panel, respectively, 
as quasi-radial modes are predominantly visible in quanti- 
ties like the density, while quadrupolar modes are prominent 
for instance in the waveform. 

The frequencies of the four main quasi-normal modes 
(the fundamental F-mode and 2 /-mode, and the harmonic 
.Hi-mode and 2 pi-mode) extracted from the migration mod- 
els and the respective values from their respective recre- 
ated post-migration equilibrium configurations excited by 
small amplitude perturbations along with the major non- 
linear self-couplings are presented in Table U for the simple 
MUn Eos and in Table [5] for the microphysical Eos. Except 
for the fundamental quasi-radial F-mode, which neverthe- 
less shows satisfactory agreement, the equilibrium models 
reproduce the frequencies of those modes to excellent ac- 
curacy, with relative differences of typically less than 5%. 
For the simple MUn EoS only in few cases the frequency 
of quasi-normal or nonlinear modes could not be identified 
(marked by dashes in Tabled, presumably because of the 



7 We note that it makes no sense to locate these post-migration 
equilibrium models in the J—fi diagram in Figs. [2] and [3] as they 
possess a small but non-negligible amount of differential rotation. 



8 For rapid rotation, the quasi-radial modes acquire a strong 
quadrupolar contribution, making an unambiguous identification 
(particularly of the first harmonic mode Hi) a demanding task. 
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Table 3. Summary of the set of stable post-migration initial models with the simple MUn EoS (family PMS) and the microphysical EoS 
(family PMM). f c ; is the central rotation frequency, A is the differential rotation length scale, and the other quantities are explained in 
Table [T] Note that as required, the quantities Mo, M and J are almost identical to the respective marginally stable models. 



Model 


M 


M 


£d 


Pc,i 


fa 


J 


A 


r p, 


iA-ci 


Ti/\Wi 




[Mq] 


[Mq] 


[10 14 g cm" 3 ] 


[10 14 g cm" 3 ] 


[kHz] 


Kc- 1 ] 








[%] 


DUO! 

r ivioi 


1 AX. 


l.OO 


1 ^ Q1 

io. yo 


lo.Oo 


u.uoy 


U.U i 


9 7 


n 
u 


QQS 
:).tc> 


n da 




i xn 

1.0U 


l.oy 


10. ol 


1 Q A A 


U.04y 


U.OU 


9 Q 
Z.£5 


u 


you 


1 dp. 


PMS3 


1.56 


1.44 


15.67 


13.33 


0.810 


0.81 


2.5 





850 


3.72 


PMS4 


1.59 


1.47 


15.53 


13.23 


0.882 


0.97 


2.7 





810 


4.77 


PMS5 


1.66 


1.54 


15.34 


13.09 


1.025 


1.27 


2.5 





740 


6.76 


PMM1 


1.35 


1.26 


16.73 


14.34 


0.498 


0.34 


2.5 





950 


1.18 


PMM2 


1.40 


1.30 


16.60 


14.23 


0.763 


0.58 


2.5 





880 


2.91 


PMM3 


1.45 


1.34 


16.51 


14.18 


0.945 


0.80 


2.5 





810 


4.73 


PMM4 


1.50 


1.39 


16.38 


14.08 


1.063 


0.99 


2.5 





750 


6.34 


PMM5 


1.54 


1.43 


16.30 


14.01 


1.118 


1.15 


2.7 





700 


7.60 


PMM6 


1.59 


1.47 


16.11 


13.87 


1.195 


1.34 


2.5 





650 


9.01 




/[kHz] /[kHz] 



Figure 10. Logarithmic frequency spectrum (i.e. Fourier transform as indicated by the hat symbol) of the time evolution of the central 
density p c (left panel) and the gravitational wave strain h (right panel) for the migration model US3 (solid line) and the associated post- 
migration equilibrium model PMS3 (dashed line, excited with an I = 0/2 trial cigenfunction for pc/h, respectively) with the simple MUn 
EoS. The vertical scale is one order of magnitude per tick mark. Several prominent quasi-normal modes and their nonlinear couplings 
are marked. The data have been scaled such that the peaks of the F-mode (left panel) and the 2 /-mode (right panel) have the same 
height for the two models US3 and PMS3. 



Table 4. Frequencies for quasi-normal modes and some of their nonlinear self-couplings for the migration models of the family US and 
the post-migration equilibrium models of the family PMS with the simple MUn EoS. /f, f2f, //jj, and f2 pi arc the frequencies of the 
fundamental quasi-radial mode, the fundamental quadrupolar mode, and their respective first harmonics, while the other modes are 
nonlinear self-couplings. For the quasi-normal modes the relative differences in frequency between the model of the US and PMS family 
are also given. For some models, particular modes could not be extracted with sufficient confidence, probably due to interaction with 
nearby other modes. 
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Table 5. Frequencies for quasi- normal modes and some of their 
nonlinear self-couplings for the migration models of the family 
UM and the post-migration equilibrium models of the family 
PMM with the microphysical MUn EoS. fp and fef are the fre- 
quencies of the fundamental quasi-radial mode and the funda- 
mental quadrupolar mode, while the other modes are nonlinear 
self-couplings. For the quasi- normal modes the relative differences 
in frequency between the model of the UM and PMM family are 
also given. For this model family, the strong damping of pulsa- 
tions prohibits a clear identification of the harmonics. 
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interaction with another sufficiently strongly excited mode 
nearby in the spectrum, while for the microphysical EoS nei- 
ther of the harmonics can be extracted unambiguously (thus 
we refrain from presenting them in Table [5]). 

A comparison of the mode spectra of the migrating 
model US3 with those of the corresponding perturbed post- 
migration equilibrium model PMS3 in Fig. \W\ also exhibits 
that the collapse associated with the dynamic migration 
in model US3 excites several nonlinear couplings between 
quasi-normal modes, with some of them being self-couplings 
like for instance the 2 • F-mode at twi ce the frequency of 
the fu ndamental quasi-radial F-mode. lAbdikamalov et alj 
(|2009l ) have also reported such strong nonlinear mode exci- 
tation in their simulations of rotating neutron star models 
that collapse to hybrid quark stars following a phase tran- 
sition to quark matter in the core. In contrast, in the only 
slightly perturbed model PMS3 these nonlinear modes are 
very weak or even not present at all. They make the same 
observation - a strong presence of nonlinear modes in the 
spectrum of the migration models and an efficient suppres- 
sion of such modes in the respective equilibrium models - 
in all members of the two model families, irrespective of 
the EoS used. We therefore speculate that in a sufficiently 
strong dynamic gravitational collapse, one can expect that 
nonlinear pulsations modes of non- negligible amplitude are 
excited, independent of the physical mechanism that acti- 
vates the collapse. 

A peculiar property of all models of family US with the 
simple MUn EoS is revealed by a comparison of the evo- 
lution of their density structure during the migration and 
the ring-down with the extracted density eigenfunction pF 
of the F-mode, which is always the mode that contains the 
largest pulsation energy (as e.g. shown by the strong peak in 
the spectrum of p c in the left panel of Fig. I10[) . We observe 
that, as expected, the density in the migrating neutron star 



changes during the evolution everywhere in the star. Sur- 
prisingly, however, the density remains practically constant 
at the two-dimensional surface where p — pi initially, as 
demonstrated for the two radial density profiles of model 
US5 along the equatorial plane and the polar axis in the left 
panel of Fig. 1111 presented at both t = and 50 ms. Appar- 
ently, while p c jumps from well below to well above p2, the 
(angular dependent) location where p — pi stays at almost 
the same radius. The plot of the corresponding density eigen- 
function profiles in the center panel of Fig. [11] then shows 
that pf has a node virtually at the same location. In the 
right panel of Fig. [TT] this behavior is summarized by plot- 
ting the maximum range during the entire evolution time of 
(rotationally symmetric) surfaces where p — p2 (closest to 
the center, with a large variation and a minimum range of 
r = as p < p2 initially in the star), where p = pi (which 
lies in a very narrow interval), and where p = p a tm (i.e. the 
boundary of the star which again changes its extent signif- 
icantly). The close proximity of the surface with p = pi to 
the nodal surface of the eigenfunction pF is striking. We 
speculate that if nodal surface is "pinned down" where the 
EoS qualitatively changes its properties at the density pi, 
which is also the location where dP/dp becomes discontin- 
uous. Apart from this attempt of explanation we are not 
able to give a compelling reason for this property of the US 
models, but it is obvious that regarding the main pulsation 
mode the migrating star decomposes into two separate do- 
mains with an interface at p = pi, each of them pulsating 
coherently, but with opposite phase. 

Similarly, we find that the models of family UM with a 
microphysical EoS exhibit a nodal surface for the F-mode 
that is approximately constant in space throughout the mi- 
gration and the subsequent ring-down, although here the 
surface expands a bit during the first phase of the migra- 
tion. However, the (again approximately constant) density 
in the neutron star model corresponding to this nodal sur- 
face now is not a distinguished location in the microphysical 
EoS, being for instance significantly smaller than pi, where 
the domain with zero pressure slope starts. Unlike in the 
models with the simple MUn EoS, in the models with the 
microphysical EoS the eigenfunction topology is apparently 
not influenced by any transition density of the EoS. 

In passing, we point ou t that in contras t to the rotat- 
ing neutron star models by iLin et all (|2006t ) that undergo 
a mini-collapse due to a phase transition in the EoS lead- 
ing to a formation of a quark matter core, our models do 
not develop significant convection after the initial contrac- 
tion. This would be easily observable both as a strong low- 
frequency contribution much below 1 kHz in the spectrum 
shown in Fig. [TO] and as a step- like pattern corresponding to 
the convection cells in the radial ro t ation velocity profiles in 
Figs. |8] and H lAbdikamalov et all (|2009l ) have shown that 
for such a phase-transiti on-induce d colla pse with the par- 
ticular EoS treatment of I Lin et al.l (|2006t ) the convection is 
caused by a significant negative entropy gradient, which is 
not present in our models. 



4.4 Emission of gravitational waves and mode 
resonance 

While the increase in density resulting from the migration 
is quite uniform for all members of each model family, the 
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r [km] r [km] x [km] 

Figure 11. Density topology of the most rapidly rotating migration model US5 of family US with the simple MUn EoS. In the left 
panel the radial profiles of the densities p e in the equatorial plane and p p along the polar axis are shown at the initial and final evolution 
time, respectively. The horizontal dotted lines mark the transition densities pi and p2, while the vertical dotted lines specify the radius 
where p c and p p equal pi apparently independent of the evolution time. In the center panel the F-mode density eigenfunction c in the 
equatorial plane and Pf,p along the polar axis is plotted. The vertical dotted lines indicate the respective nodes. In the right panel we 
show the range (symbolized by arrows) for the radial location of the surface where the density p equals p2, pi, and the atmosphere cutoff 
density p a tm during the evolution, as well as the nodal surface of pp. Here x = r s'm8 and z = r cos 9 are the Cartesian coordinates. 



strength of gravitational wave emission, which is sensitive to 
the rotation frequency, is expected to vary as the initial value 
of / and the change A/ during the collapse varies among 
different models. The dependence of the waveform on the 
rotation rate is clearly revealed in the gravitational radia- 
tion waveform plots presented in Figs. [T2l and 1131 for three 
selected models of each of the two families US and UM, re- 
spectively. The according (absolute) peak amplitudes |/i| max 
of the waveforms are summarized in Table[2] The main emis- 
sion modes are the 2 /-mode and, for rapidly rotating models, 
also the F-mode. 

It is remarkable for this migration scenario that |/i| max 
apparently does not monotonically depend on either /, or 
Ti/|Wi| (see e.g. the plot in the top panels of Fig. 1141 for 
model family US and Fig. 1 151 for model family UM). An in- 
structive example is also the relatively slowly rotating model 
US1 (see left panel of Fig. I12[) , where the maximum wave- 
form amplitude during the initial collapse phase of the mi- 
gration, which lasts at most 2 ms, is around 10~ 22 at a 
distance of 10 kpc from the source. Considerably later, at 
t ~ 5 ms, already well in the ring-down phase, the amplitude 
gradually increases until it saturates at times after ~ 10 ms 
with a maximum |/i| max ~ 5 x 10 -22 . Finally, after that the 
amplitude of the quasi-periodic waveform's envelope decays 
as expected (see also the qualitatively similar waveform of 
model UM1 in the left panel of Fig.[T3J. 

A very similar behavior was also found in recent simu- 
lations of rotating neutron st ar models which collapse due 
to an EoS phase transition (|Abdikamalov et ail 12009 ) . In 
that work it was shown that nonlinear modes can vigorously 
boost the emission of gravitational waves if they are in res- 
onance with an efficiently emitting quadrupolar mode like 
the 2 /-mode, provided that both the frequencies of those two 
modes are in close proximity and the nonlinear mode con- 
tain s a substantia l amount of pulsation kinetic energy (se e 
also lEardlevI [19831 ; iDimmelmeier, Stergioulas fc Fontll2006l) . 



While in the models of lAbdikamalov et al. I (|2009h the non- 
linear mode which drives the mode resonance is the 2 • F- 
mode, the mode to be considered here is the 3 • F-mode for 



the models with the simple MUn EoS and the 2 • F-mode for 
the ones with the microphysical EoS, as suggested by the 
frequencies given in Tables U and [S] respectively. 

Indeed Figs. Q3] and \T5\ show the connection between 
the proximity of these potentially resonating modes and the 
maximum waveform amplitude |/i| max for all migrating mod- 
els considered here. By evaluating |/i| max on the one hand 
only in a time window close to the time of maximum com- 
pression (i.e. at the end of the dynamic initial collapse phase 
before mode resonance can develop; dashed line) and on the 
other hand for the full waveform (solid line), the enhance- 
ment of gravitational wave emission by resonance can be 
quantified. For the models with a microphysical EoS of fam- 
ily UM in Fig. 1151 both the absolute and the relative am- 
plification of the waveform amplitude (see top and center 
panel, respectively) is clearly strongest when the 2 ■ F-mode 
is closest in frequency to the 2 /-mode. This happens at low 
rotation, while for rapid rotation we detect no imprint of 
mode resonance in the waveform and |/i| ma x obtained close 
to the time of maximum compression is the global maximum. 
The models with a simple MUn EoS presented in Fig. [TJ] re- 
veal a slightly more complicated behavior. Again the slowest 
rotator exhibits the strongest relative increase of |ft| ma x by 
mode resonance (see center panel), albeit from a very small 
absolute level (see top panel). As rotation grows and the 
frequencies of the 2 • F-mode and the 3 • F-mode approach 
each other, the effects of the resonance on |/i| max become ap- 
parent. However, in the rapid rotation limit mode resonance 
again ceases (as for the models of family UM), although 
in family US the frequencies of the two relevant modes are 
now in close proximity. Here the waveform again possess a 
regularly damped sinusoidal pattern. 

For such rapidly rotating models, we find that the peaks 
associated to the 2 /-mode and the 3 • F-mode (for family 
US) or the 2 • F-mode (for family UM) are clearly sepa- 
rated also during the interval from 5 to 15 ms. Additionally, 
the respective nonlinear self-coupling of the F-mode is ex- 
cited to a much lesser degree than in slower rotating mod- 
els with mode resonance, and thus carries less transferable 
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Figure 12. Time evolution of the gravitational wave strain h for three migration models of family US with the simple MUn EoS. 




t [ms] 

Figure 13. Time evolution of the gravitational wave strain h for three migration models of family UM with the microphysical EoS. 



energy. Furthermore, rapid rotation also causes the regular 
fundamental quasi-radial F-mode to acquire such a strong 
quadrupolar nature that it becomes the main gravitational 
wave emitting mode, surpassing the 2 /-mode in peak height 
in the waveform spectrum (see also the discussion in Sec- 
tion I4.5|l . Therefore, even if for such models the 2 /-mode 
was somewhat amplified, this would not be identifiable in 
the spectrum, which is dominated by the F-mode. We also 
note that the inefficacy of this t ype of mode resonance for 
rapid rotation was also found bv lAbdikamalov et al.l (|2009l ) 
in their simulations of the collapse of neutron star models 
to hybrid quark stars. 

We also note that although in the dynamic collapse 
phase of a migration the self-couplings of the F-mode, which 
are responsible for transferring energy into the emitting 2 f- 
mode, are much more violently excited, they are also present 
in perturbed equilibrium models. Consequently, we observe 
the effects of mode resonance not only in some of the mi- 
grating models of family US and UM, but also in the cor- 
responding equilibrium models of family PMS and PMM, 
if these are subjected to a small amplitude initial pertur- 
bation according to Eq. (|19|) . which favors the excitation 
of quasi-radial / = modes. As expected, the quadrupolar 
I = 2 perturbation of Eq. (|20|l excites self-couplings of the 



F-mode only weakly, and thus in that case mode resonance 
does not occur in any of the equilibrium models. 



4.5 Detectability prospects of gravitational waves 

Due to the long quasi-periodic emission time the gravita- 
tional waves emitted by a migration event is a promising 
source for gravitational wave detectors of both interferome- 
ter or resonant type. Assuming that the strong post-collapse 
oscillations are not damped by any other physical mecha- 
nism than dissipation of kinetic energy by shocks, the damp- 
ing we observe for the models of family US suggests that the 
effective total emission time for gravitational waves can be 
much longer than the evolution time of our models (which 
is tf = 50 ms), i.e. last for hundreds of pulsation cycles. In 
contrast, the models of family UM are so strongly damped 
that the emission of gravitational waves has essentially died 
out at the end of the evolution time at t — 50 ms. 

To assess the prospects for detection by current and 
planned interferometric detectors, we calculate character- 
istic quantitie s for the gravitational wave signal following 
iThornel (| 1987T ) . Applying a Fourier transform to the dimen- 
sionless gravitational wave strain h, we obtain 
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Table 6. Detection prospects for gravitational waves. f c is the characteristic frequency, h c is the integrated characteristic gravitational 
wave signal strain and SNR is the signal-to-noise ratio, each given for the LIGO detector, the VIRGO detector, and the Advanced LIGO 
detector in broadband mode. All quantities assume a total emission time of tf = 50 ms and are dependent on the rms strain noise /i rms 
of the detector. 



Model 


/c.LIGO 


/c.VIRGO 


/c.Adv. LIGO 


he, LIGO 


frc.VIRGO 


he, Adv. LIGO SNRliGO 


SNRvirgo 


SNRAdv. LIGO 




[khz] 


[khz] 


[khz] 


" lO" 21 




" lO" 21 




" lO" 21 ' 


[at 10 kpc] 


[at 10 kpc] 


[at 10 kpc] 










at 10 kpc 




at 10 kpc 




at 10 kpc 








US1 


1.768 


1.786 


1.798 


2.6 


2.6 


2.6 


0.3 


0.5 


9.2 


US2 


1.225 


1.321 


1.432 


10.4 


11.0 


11.7 


2.1 


3.4 


54.9 


US3 


1.574 


1.651 


1.728 


31.6 


32.9 


34.1 


4.4 


7.6 


127.6 


US4 


0.924 


1.006 


1.129 


9.3 


9.8 


10.7 


2.9 


1.1 


65.0 


US5 


0.718 


0.766 


0.834 


5.8 


6.0 


6.4 


2.6 


3.7 


51.8 


UM1 


2.024 


2.031 


2.041 


18.6 


18.7 


18.7 


1.8 


3.2 


56.6 


UM2 


1.928 


1.952 


1.983 


22.0 


22.2 


22.5 


2.3 


1.0 


70.6 


UM3 


1.236 


1.297 


1.384 


7.0 


7.2 


7.6 


1.4 


2.3 


37.0 


UM4 


1.130 


1.190 


1.282 


7.9 


8.2 


8.6 


1.8 


2.9 


45.7 


UM5 


1.047 


1.110 


1.202 


7.5 


7.8 


8.3 


1.9 


3.1 


47.3 


UM6 


1.007 


1.078 


1.178 


8.8 


9.2 


9.8 


2.4 


3.8 


57.3 



10 



"i 1 r 




i 1 r 




simple MUn EoS 

4 ' 5 ' 6 




2.2 
2.1 
2.0 
1.9 
1.8 
1.7 
1.6 



_ 1 1 1 1 1 1 1 1 1 1 


i 1 i 1 i 1 i 1 i _ 


- %- »-^ 










\ 


- microphysical EoS 

",1,1,1,1,1 


\ 

\ 

• — 

,1,1,1,1,' 







T./\W. \ [%] 



4 5 6 
T./\W.\ [%] 



9 10 



Figure 14. Dependence of the maximum gravitational wave 
strain |/i| max (top panel), the relative amplification of the maxi- 
mum gravitational wave strain A | h | max ro i (center panel) , and the 
frequency / for potentially resonating modes (bottom panel) on 
the initial rotation rate 7] / 1 Wi | for the migration models of fam- 
ily US with the simple MUn EoS. |/i| max is given with (solid line) 
and without (dashed line) taking into account the effects of mode 
resonance, and A|/i| max rc i is the relative increase in |/i| max due to 
mode resonance. The frequencies of the 2 /-mode and 3 • F-mode 
are marked with a solid and dashed line, respectively. 



Figure 15. Dependence of the maximum gravitational wave 
strain |/i| max (top panel), the relative amplification of the maxi- 
mum gravitational wave strain A|/i| max ro i (center panel), and the 
frequency / for potentially resonating modes (bottom panel) on 
the initial rotation rate T\ / 1 Wi for the migration models of family 
UM with the microphysical EoS. |/i| max is given with (solid line) 
and without (dashed line) taking into account the effects of mode 
resonance, and A|/i| max rc i is the relative increase in |/i| max due to 
mode resonance. The frequencies of the 2 /-mode and 2 • F-mode 
are marked with a solid and dashed line, respectively. 
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h= / e 27V ' ft hdt. (21) 

J — oo 

One can compute the (detector dependent) integrated char- 
acteristic frequency as 

*-(f£'«)(jHH' 

and the dimensionless integrated characteristic strain as 

fc>= (3 y ^{h 2 )fdf) : (23) 

where is the power spectral density of the detector and 
She = Sh(fc)- We approximate the average (h 2 ) over ran- 
domly distributed angles by h 2 , assuming optimal orien- 
tation of the detector. From Eqs. (|22l I23[) the signal-to- 
noise ratio can be calculated as SNR = h c /[h Ilna (f c )], where 
/i rms = \fJB~h is the value of the rms strain noise for the 
detector (which gives the theoretical sensitivity window). In 
Table[6]we summarize the values of / c , h c and the signal-to- 
noise ratio for all the migrating models of the two families US 
and UM for the currently operating LIGO and VIRGO de- 
tectors and for the future Advanced LIGO detector in broad 
band operation mod^S Here we consider a source inside our 
own Galaxy at a reference distance of 10 kpc. 

The fraction of neutron stars that undergo a migration 
process at some stage in their lifetimes is rather uncertain. 
In this work we restrict ourselves to considering isolated neu- 
tron stars. Depending on the selected gravitational mass and 
the EoS, such a neutron star must have spun down to some- 
where between 50 and 800 Hz to be subject to a migration 
(see Table[T] and Figs. [2] and [3}, provided it has an EoS that 
is similar to one of our choice and shows the back bending 
phenomenon. 

Apparently, isolated millisecond pulsars with low mag- 
netic field, which have been recycled by accretion, are no 
appropriate candidates, as they spin down too slowly, on 
a Hubble time scale. However, there exists a another class 
of neutron stars, which are thought to rotate initially at a 
millisecond period, and those are magnetars. In order to gen- 
erate the required surface magnetic field of roughly 10 14 to 
10 15 G, the proto-neutron star born in a type II supernova 
explosion must rotate at ~ 500 Hz. Due to strong magnetic 
breaking, it slows down to a few Hz after only about 10 3 y. 

The most dramatic spin-down of a magnetar occurs dur- 
ing the first years of its life, which is the period when the mi- 
gration phenomenon could occur. So far ten magnetars were 
detected in Galaxy: six anomalous X-ray pulsars and four 
soft-gamma repeaters. Magnetars are active and therefore 
observable for about 10 3 - 10 4 y. This means that their birth 
rate is not much smaller than that of "ordinary" radio pul- 
sars (with a lifetime of ~ 10 6 y) , which in turn is comparable 
to the type II supernova rate in a big spiral galaxy like ours. 
Consequently, we may expect about one migration event in 
our Galaxy per century, which is a very low rate. Consistent 

9 Note that both in Table[6]and in Fig. ll6l thc gravitational wave 
characteristics have been evaluated for a total emission time of 
t{ = 50 ms, where the waveform amplitude of at least the models 
of family US has not yet decayed to zero. Therefore, for these 
models the values for h c and the signal-to-noise ratio shown here 
are actually a lower bound. 



with this simple argument, an analysis based on observa- 
tional data yields a very uncertain birth rate betwee n 0.2 
and 6 per century (|Muno et al.ll2007l ; lGill fc Hevlll200 7T ) . The 
magnetar birth rate rises to about one p er year if the whole 
Virgo cluster of galaxies is considered (|DaH'Osso fc Stellal 
l2007n . 

Still, if such a migration event did occur in our Galaxy, 
for current interferometric detectors of the LIGO class and 
assuming an emission time if = 50 ms, all of our models 
except the slowly rotating one US1 have a signal-to-noise 
ratio very close to or above 1. Model US3, which is both 
moderately rapidly rotating and whose gravitational wave 
emission is amplified by mode resonance (see center panel 
of Fig. I12|l . features the largest signal-to-noise ratio of all 
investigated models at about 4 for current LIGO and 8 
for current VIRGO. For the Advanced LIGO detector, the 
signal-to-noise ratio lies comfortably above 10 for all mod- 
els (again except model US1 which still yields SNR = 9) 
and reaches ~ 130 in model US3. It is remarkable that h c 
(and correspondingly, also the signal-to-noise ratio) for the 
model family UM with the microphysical EoS is bounded 
by a rather narrow interval between about 7 and 22 x 10 -21 
at 10 kpc, although the initial rotation rates of the models 
span a wide range (see Table [TJ . 

For substantially increasing the event rate, it would be 
necessary for the detector to be sensitive to signals com- 
ing from distances out to the Virgo cluster, at ~ 20 Mpc. 
However, at this distance the signal-to-noise ratios for our 
models given in Table [6] must to be scaled down by a fac- 
tor of 2000 and consequently drop to well below 1 even 
for Advanced LIGO. Therefore, as for gravitational wave 
signals from supernova c ore collapse (see the discussion in 
iDimmelmeier et alj|2007h . the second generation detectors 
will improve the signal-to-noise ratio of a local event alone, 
but will not increase the event rate much on account of 
the inhomogeneous galaxy distribution in the local region 
of the Universe. Only third generation detectors will have 
the required sensitivity in the kHz range to achieve a robust 
signal-to-noise ratio at the distance of the Virgo cluster. 

Note that for most of the models the integrated char- 
acteristic frequency f c given in Table [6] is not very close to 
either of the two main gravitational wave emission frequen- 
cies /f and f2f. This is because f c reflects the frequency 
dependence of the sensitivity of the detector, because non- 
linear mode couplings and higher-order linear modes also 
contribute to the gravitational wave signal (although with 
lower relative amplitudes) and, most importantly, because 
for many models the gravitational wave power spectrum of 
the signal exhibits two strong peaks from the F-mode and 
the 2 /-mode. Based on the peak height in the spectrum of 
the gravitational radiation waveform, we find that the fun- 
damental quadrupolar 2 /-mode is the strongest emitter for 
the slowly and moderately rapidly rotating models, followed 
by the fundamental quasi-radial F-mode. For rapid rotation, 
this order of importance is reversed. 

The detector dependence of f c and h c is also illustrated 
in Fig. 1161 where the locations of the gravitational wave sig- 
nals for all models are plotted relative to the rms strain noise 
/i rms of the current LIGO and VIRGO detectors (left and 
center panels, respectively) and the future Advanced LIGO 
detector (right panel), all for a distance to the source of 
10 kpc. Model US3, which both has significant rotation and 
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Figure 16. Location of the gravitational wave signals for the migration models from family US with the simple MUn EoS (filled circles) 
and family UM with the microphysical EoS (open circles) in the h c —fc plane relative to the sensitivity curves (i.e. the rms strain noise 
forms) of the LIGO detector (left panel), VIRGO detector (center panel), and Advanced LIGO detector in broadband mode (right panel), 
all at a distance of 10 kpc. Models whose gravitational wave emission is amplified by mode resonance effects are marked with a cross. 



exhibits mode resonance (models where more resonance am- 
plifies gravitational radiation emission are marked by crosses 
in Fig. 116ft . clearly sticks out with the largest values of h c . 
In general, models with initially moderate rotation, whose 
waveform amplitudes are in addition boosted by mode reso- 
nance effects, are found in the upper right quarter of the di- 
agram, i.e. they feature both a high characteristic frequency 
and a large characteristic strain. Model US1, which is by far 
the slowest rotator, shows the smallest value for h c , albeit 
at a comparably high f c . 

From Fig.[l6]it is clear that a possible detection is ham- 
pered by the adverse characteristics of the detector sensitiv- 
ity curve in the frequency interval relevant for our models, 
where the strain noise h lms rises steeply with frequency. Pro- 
posed future detectors with a flat h TUls curve up to frequen- 
cies in the kHz range could eliminate this problem. 

Interestingly, if a larger selection of models that are 
evenly distributed in the parameter space of initial rotation 
rate were simulated, their distribution in the h c —f c plane 
would m ost probably very well ap proximate the pattern 
found by iDimmelmeier et al.l l)2008l ) for rotating models of 
collapsing stellar cores in a supernova event. In order to fil- 
ter out the effects of mode resonance, which for some models 
amplifies the waveform amplitude in the ring-down following 
the collapse phase and which cannot be simply subtracted 
from the time integrated characteristic strain h c , in Fig. 1171 
instead of h c we plot the maximum gravitational wave strain 
|/i|max around the time of maximum compression (i.e. the 
most dynamic moment of the collapse phase) before any 
possible waveform amplification by mode resonance. 

As for he in Fig. 17 of IDimmelmeier et al.l <|200Sh . 
starting from very slowly rotating models |/i| max and (not 
as strongly) f c increase as rotation grows (see arrow 1 
in their figure and our Fig. I17[l . then only the wave- 
form amplitude rises at practically constant f c (arrow 2 
again in both figures). We recall that for these mod- 
els the 2 /-mode is the strongest emitter, which depends 
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Figure 17. Dependence of the maximum gravitational wave 
strain /i| max on the characteristic frequency f c for the migration 
models of family US with the simple MUn EoS (unfilled circles) 
and family UM with the microphysical EoS (filled circles). If mode 
resonance influences gravitational wave emission (models denoted 
with a cross), |/i|max also is computed around the time of maxi- 
mum compression before the amplification of h takes effect. The 
decrease of |/i| max without resonance effects is marked by dot- 
ted lines. The arrows indicate the qualitative effect of increasing 
initial rotation rate Ti/|Wi|. 



only weakly on rotation (see Tables [4] and [5] and also 
IDimmelmeier. Stergioulas fc F ont 2006), which explains the 
approximate constancy of f c in this regime. This weak im- 
print of rotation is in stark contrast to e.g. a rigidly rotating 
bar, where the emission frequency increased linearly with 
rotation frequency. 

Finally, for rapidly rotating models which run into 
the centrifugal barrier, |/i| max saturates at a constant or 
even slightly declining level with significantly decreasing fre- 
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quency f c (arrow 3 again in both figures). What we do not 
observe here is the rather separate group of extremely rotat- 
ing models, whose evolution and waveforms are dominated 
by centrifugal forces (marked by 4 in their figure). In addi- 
tion, in our case mode resonance effects would bulge out the 
"knee" for moderately rotating models in the upper right 
corner of the curve towards higher values of h c (or in our 
case |/i| max ), as seen in Fig. 1171 where for such models the 
maximum amplitude for the entire waveform (marked by 
crosses) is higher than if only taken in the period around 
maximum compression. 



5 SUMMARY AND CONCLUSIONS 

We have studied the mini-collapse of rotating neutron star 
models induced by a dynamic back bending instability. This 
instability is caused by the softening of the equation of state 
resulting from a phase transition in the dense core. The col- 
lapse leads to a migration of the neutron star from an unsta- 
ble equilibrium configuration to a stable one after a series of 
damped pulsations. The migration starts when the configu- 
ration reaches the local minimum of the total angular mo- 
mentum in the back bending diagram. We have investigated 
two families of models, one with the simple analytic MUn 
equation of state and one with a tabulated microphysical 
equation of state that includes effects of a phase transition 
due to kaon condensation. 

If the rotation profile of the neutron star remained uni- 
form during the collapse, the neutron star would simply mi- 
grate horizontally (i.e. at constant total angular momentum) 
from the unstable branch to the stable branch (which are 
separated by the unstable segment) in the angular momen- 
tum vs. rotation frequency diagram. However, our simula- 
tions show that the neutron star actually develops a differ- 
ential rotation profile. Still, we find that the final degree of 
differential rotation is not large and that the new equilib- 
rium configuration is similar to the corresponding one on 
the stable branch in the back bending diagram under the 
assumption of rigid rotation. 

The dynamic collapse in the migration process creates 
a whole spectrum of coupled pulsation modes, which we ex- 
tract directly from the numerical results of our fully nonlin- 
ear simulations. In addition, we have also constructed equi- 
librium models which closely resemble the post-migration 
configurations, and evolved them after applying a small 
perturbation. Using both approaches, we have identified 
the characteristics of the main modes excited by the mini- 
collapse, and find evidence for resonance between the fun- 
damental quadrupolar mode and nonlinear self-couplings of 
the fundamental quasi-radial mode. Our results show that 
the amplitude and waveform of gravitational radiation de- 
pends strongly on the mass of the star and the equation of 
state used. Namely, models with the microphysical of state, 
which features a density jump associated with a first-order 
phase transition, suffer from much stronger damping than 
models with the mixed-phase MUn equation of state. 

In the case of a detection with current or future gravita- 
tional wave interferometer detectors (like VIRGO or LIGO), 
such imprints of the model parameters on the waveform 
properties could be exploited to observationally constrain 
the unknown equation of state for dense matter. For a Galac- 



tic event, whose frequency of occurrence should be compa- 
rable with the local type II supernova rate, the predicted 
gravitational wave signal from such a dynamic migration lies 
above the designed sensitivity curves of current km-size de- 
tectors. With future, more sensitive detectors, signals from 
events within the Virgo cluster of galaxies could also come 
into range. Then the corresponding event rate could pos- 
sibly increase to about one per year. In addition, we find 
that in some cases the gravitational wave emission is con- 
siderably enhanced due to the said mode resonance, which 
transfers energy from the strongly excited but weakly radi- 
ating quasi-radial mode to the initially weakly excited but 
strongly radiating quadrupolar mode. 

Particularly promising candidates for a dynamic migra- 
tion are young magnetars, as they have an initial rotation 
rate in the kHz range and spin down very rapidly. While 
the prospective astrophysical phenomenon we study in this 
work seems to be rare, it has nevertheless a potential of pro- 
ducing a very characteristic signature in the form of a pe- 
culiar burst of gravitational radiation. The damping of the 
post-migration pulsations, which is directly reflected in the 
gravitational wave signal, is even sensitive to the character 
of the phase transition softening, i.e. whether it occurs with 
or without a jump in density. Therefore, if chance allows 
us to detect such a dynamic migration event in our Galaxy, 
this will possibly be a source of rich information about the 
structure of the superdense neutron star core. 

We point out that in our simulations the effects of the 
magnetic field, which for magnetars has a typical strength 
of B ~ 10 14 -10 15 G, are neglected. In particular the impact 
on the stellar shape is discarded. However, the typical pres- 
sure in neutron star cores, where the phase transition in the 
equation of state that causes the back bending phenomenon 
occurs, is 10 34 -10 35 dyn cm -2 (see Figure [TJ, while with 
B 2 /{8n) ~ 4 x 10 29 (B/10 15 G) 2 the stress due to the mag- 
netic field is many orders of magnitude lower. Therefore, the 
influence of even a strong magnetic field on the core struc- 
ture is small, and consequently we expect only secondary 
effects on the back bending instability. 

We close by pointing out that the results presented here 
should be considered as a first step towards a better under- 
standing of the neutron star interior and dynamical phe- 
nomena like migration from unstable to stable configura- 
tions. Due to the presently very limited knowledge about the 
equation of state for dense matter, some of the aspects of the 
physical processes underlying the investigated scenario are 
necessarily simplified. The treatment of the equation of state 
as a representation of matter in the mixed phase experienc- 
ing a phase transition is rather crude, especially when one 
considers the behavior of fluid elements undergoing succes- 
sive stages of compression and decompression, and changing 
the proportions of the phases. In particular, our modeling 
neglects any possible effects from local heating, or from the 
creation and subsequent emission of neutrinos during the mi- 
gration collapse and the subsequent bounces. On the other 
hand we have shown that the time scale of the mini-collapse 
is short enough to neglect non-equilibrium and bulk viscos- 
ity effects. 
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APPENDIX A: NON-EQUILIBRIUM EFFECTS 
IN THE EQUATIONS OF STATE 

In its quasi-stationary evolution phase, the rotating neutron 
star slowly proceeds along the stable branch. Here the slow 
angular momentum loss results in an according compression 
of the neutron star matter. When it enters the dynamic mi- 
gration period, the phase transition in the stellar core leads 
to a much more rapid initial contraction, followed by several 
re-expansion and compression cycles (see Figs. [6] and [7}. In 
both regimes the composition of matter has to re-adjust to 
the change in density and pressure to keep matter in thermo- 
dynamic equilibrium. This process of equilibration of com- 
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position proceeds via nuclear reactions, of which the ones 
involving weak interaction (changing lepton number and/or 
strangeness per baryon) are the slowest. We now investigate 
these relaxation processes and their impact on the EoS of 
compressed matter separately for the two regimes. 

When the neutron star evolves on the stable branch 
(marked by solid lines in Figs. [2] and [3}, it move slowly 
downwards along the branch on a characteristic secular time 
scale t scc = pc/pc oc (J/|J|) _1 . This time scale is so long 
that equilibration nearly catches up with compression, and 
any deviations from equilibrium are tiny. Consequently, the 
EoS of matter can be well approximated by assuming full 
equilibrium. 

As soon as the marginally stably equilibrium config- 
uration is reached, further loss of J triggers the dynamic 
migration, associated with a compression of matter on the 
dynamic time scale Td yn < 1 ms. Now matter moves off weak 
equilibrium and the actual EoS stiffens compared to the fully 
equil i brated one (see e.g. iGourgoulhon. Haensel fc Gondekl 



1 19951 ; lHaensel. Levenfish fc Yakovlevl |20o3 T and references 
therein). However, a deviation from equilibrium accelerates 
(in a strongly nonlinear manner) the rate of equilibration 
processes, which tend to return the matter towards equilib- 
rium following le Chatelier's principle. Consequently, also in 
the migration phase (i.e. the rapid initial contraction and the 
subsequent pulsations) the EoS is quite close to the equilib- 
rium one. Strictly speaking, the deviation from equilibrium 
is such that as to make the rate of equilibration equal to 
the time scale of dynamic compression. Still, the effect on 
the dynamics is small, and is thus being neglected in our 
calculation. 



APPENDIX B: PULSATION DAMPING IN THE 
MICROPHYSICAL EQUATION OF STATE 

The damping of the post-migration pulsations observed in 
the models of family UM with the microphysical EoS is much 
stronger than in family US with the simple MUn EoS (cf. 
Figs.[7]and[5] respectively). It is caused by a particular, very 
efficient damping mechanism that occurs if matter moves 
across a first-order phase trans ition in such a type of EoS 
due to quasi-radial pulsations l|Bisnovatvi-Kogan fc Seidovl 
1 1984 lHaensel. Zdunik fc Schaefferlll989T ) In the following, 
we explain that mechanism and the relevant energy and time 
scales in our models in more detail. For simplicity, without 
loss of generality we do not take into account effects of ro- 
tation and only consider spherically symmetric Newtonian 
models. 



Bl Radial pulsations of the dense core 

First we consider a rather simplified model of a pulsating 
compact star with an EoS that exhibits a phase transition 
and leads to a dense core in the center of the star. We ne- 
glect finite compressibility of matter outside the core (except 
for the phase transition), so that the star's density there is 
simply p\. 

Conservation of mass then leads to a relation between 
the pulsation amplitude Sr c of the dense core radius and the 
radial velocity u r >r c of the pulsations outside the core, 
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Figure Bl. Time evolution of the radius r c , P = P1 of the dense 
core in the equatorial plane for the slowest rotating model UM1 
of family UM with the microphysical EoS. The pulsation period 
of r e ,p= P1 is roughly t p = 1 ms. The pulsation amplitude 8r c is 
damped to 1/5 of its initial value of about <5r c o = 0.55 km in 
approximately rj/g = 6 ms. The equatorial radius of the dense 
core finally settles at the equilibrium value r c = 4 km. 
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where r c is the unperturbed equilibrium radius of the dense 
core, r is the unperturbed radial coordinate, u> is the angular 
pulsation frequency, and A = pa /pi quantifies the jump in 
density at the core edge due to the phase transition in the 
EoS. The phases of 8r c (t) and v r>Tc (t) differ by tt. As the 
core matter is incompressible, the velocity v r<Tc at radii 
smaller than r c vanishes. 

Under our simple model assumptions, the average ki- 
netic energy contained in the radial pulsations is then given 
by 



E v = 7T / drr 2 pi vj5 ; 



= ^P!r c 3 o; 2 (A-l) 2 (l-^)(«Sr c ) 2 , 



(B2) 



where r s is the unperturbed equilibrium radius of the star. 

The stable equilibrium initial model PMM1 that cor- 
responds to the post-migration state of model UM1 has a 
stellar radius r Sj o ~ 12 km, while from Fig. IB1I we can ex- 
tract typical values <5r Ci o ~ 0.55 km and r Ci o ~ 4 km at the 
start of the post-migration pulsationf°l. For our microphys- 
ical EoS the density discontinuity corresponds to A = 1.62 
(see Section niTI and Fig. [Q). Thus, from Eq. (1B2|) we obtain 
for the initial pulsation energy: 



E 



p,0 



3.7 x 10 48 erg. 



(B3) 



10 For the other models of family UM we obtain similar values for 
<5r c ,o and r Cj o and r a ,o- In Fig. IBll we show the evolution of r p — pi 
only in the equatorial plane, as both the qualitative behavior 
and the typical scale of this quantity does not depend much on 
latitude. 
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B2 Dissipation from a first-order phase transition 
and associated damping time scale 

We now consider the situation where the matter on 
both sides of r c is fully equilibrated, i.e. realizes a 
minimum of the appropriate t hermodynamic potential 
ijBisnovatvi-Kogan fc Seidovll 19841 ). The kinetic energy asso- 
ciated to the pulsations is dissipated during the first-order 
phase transition, as matter passes through the boundary 
of the dense corj 1 " 1 ! Actually, the dissipation occurs during 
the compression phase, which is the first half of the pulsa- 
tion period t p , and takes place within the shock front that 
forms at the dense core edge. In the thin layer between the 
two phases, the sound speed c s « [(P2 — Pi)/{£2 — £i)] 1//2 
is rather low, because the large density jump is associated 
with only a very small pressure increase. Therefore the mat- 
ter flow within this layer is supersonic. In reality, the dissi- 
pation within the shock heats up the matter that flows into 
the dense core. However, in our approximation we neglect 
thermal effects in the EoS, and thus dissipation manifests 
itself only via the loss of kinetic energy in the matter flovF^I. 

Th e average energy dissipated during one oscillation pe- 
riod is (|Bisnovatvi-Kogan fc Seidovll 1984 ) 



and using Eq. (|B5|I . Eq. (jBTj can be resolved for 5r c : 



(B4) 



This dissipation results in a damping of pulsations and the 
decrease of Sr c . 

Assuming that all terms in Eq. (|B4|l except 8r c are con- 
stant, we obtain for the average dissipated power 

3 



E d 



z Ed,t p /t p =Ed 



5r c 
<5r c ,0 



(B5) 



Evaluating Eq. (|B4[) at the initial pulsation time and using 
the approximate value t p w 1 ms extracted from Fig. IB1I 
together with the other estimates from Appendix IBll and pi 
for the microphysical EoS, the power dissipated during the 
first pulsation period is 



E d , « 3 x 10 48 erg ms" 1 . 



(B6) 



We now demand that the change in pulsation energy 
is entirely caused by the above dissipation mechanism, and 
thus we arrive at 



dE 
~~di 



- = -Ed. 



(B7) 



Rewriting Eq. (|B2|) for the average kinetic energy associated 
to the pulsations as 



E p =E Pi o 



Sr c 



(B8) 



11 Unlike in the estimate in Appendix IB II here the density out- 
side the core need not be constant. 

12 As explained in Footnote [4] in our simulations we introduce 
a small pressure gradient between the two phases, and thus the 
flow there remains subsonic. Nevertheless, the largo density jump 
leads to numerical dissipation of comparable strength. We have 
also performed tests with almost zero pressure gradient and thus 
supersonic flow across the density jump, and find similar damp- 
ing time scales of at most a factor 2 shorter. However, for numer- 
ical reasons here strong high-frequency noise superimposes the 
damped pulsations. 



d{8r c f 



1 Ea,o 



(Sr c ) 3 dt 5r c>0 E P}0 ' 

This can be straightforwardly integrated to 



(B9) 



Sr c = Sr c ,o 1 



-Bd.U t 

2E 



p,0 



Sr c 



2E 



p,0 



Ed.o t 



(BIO) 



With this we have an explicit expression for the time 
evolution of the pulsation amplitude 5r c of the dense core, 
which we can apply to get an estimate for the damping time 
scale of the pulsations. Choosing for instance t in Eq. f|B10f) 
as the time Ti/5 when the radial perturbations Sr c of the 
dense core have dropped to <5r Cj o/5 and using the values 

from Eqs. (|B3IIB6|I for_ZJ Pi o and_Bd,o, respectively, we obtain 
r i/5 ~ 12 nis. The value T1/5 w 6 ms read off from Fig. IB1I 
is shorter, but of the same order of magnitude (~ 10 ms) as 
this crude analytic estimate. 

Although this treatment is somewhat simplistic, it 
demonstrates that the compression of matter passing 
through the phase transition in the microphysical EoS can 
very efficiently dissipate pulsation energy on a time scale 
that is consistent with our observations for the family UM 
of migration models. In contrast, as apparent from Fig. [5] 
the typical damping times of the models of family US with 
the simple MUn EoS are much longer, as there the above 
efficient dissipation mechanism is not at work. 



B3 Non-equilibrium effects and bulk viscosity 
damping 

We now relax the condition of complete equilibration. The 
phase transition at the core radius r c must then be treated 
by taking into account the slowness of weak interaction 
processes. Consequently, kaon condensation will proceed 
off thermodynamic equilibrium. However, the rate of non- 
equilibrium reactions increases very steeply with the devi- 
ation from equilibrium. Therefore, also in this case E d is 
not significantly smaller than the upper bound obtained in 
Appendix IB2I under the assumption of full equilibration. 

In addition, non-equilibrium processes produce entropy 
in the form of heating up the matter. This effect can be 
represented by bulk viscosity, corresponding to an additional 
average dissipation power Z?d, tw- 
in the case of a mixed quark- nucleon core, the effi- 
cient damping mechanism described in Section |B2I does not 
work, as there exists no single phase interface with a large 
jump in density. Then the bulk viscosity generated by non- 
equilibrium, strangeness changing reactions in the quark 
component, 

u + d^s + u, (Bll) 

becomes the main damping mechanism fsee lMadsenlll99"9l . 
and references therein). 

If one considers the mini-collapse of an unstable neutron 
star caused for instance by dynamic migration as in this 
work rather than a supernova core collapse, the relative post- 
collapse pulsation amplitude |<5r s |/r s is typically less than 
6 x 10~ 2 T/10 9 K. Therefore, bulk viscosity does not depend 
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on Sr s and the associated damping exhibits the usual form 
(|Madsenlll999l ) 

5r s = <5r s , e- t/Tbv <-, (B12) 

where Tb v , q is the damping time scale. For a strange quark 
mass of 100 MeV, a Fermi energy of the down quark of 
300 MeV, and choosing u p = 2-k /t v « 2 x 10 4 s -1 on e gets, 
after transforming a formula given in <|Madsenlll999h . 

r bv , q « 0.1 (T 9 )- 2 B. (B13) 

If the core, however, consists of a kaon condensate, then 
the bulk viscosity is generated by the non-equilibrium pro- 
cesses 

nUp + R- (B14) 

that also change strangeness. For non-superfluid nucleons, 
this bulk viscosity is independent of temperature and causes 
a damping of the fundamental radial mode pul sations on a 
time scale (|Chatteriee fc Bandvopadhvay||200'it ) 

rb ^ 100 (ofe) ' s ' (B15) 

which is much longer than the corresponding bulk viscosity 
damping time scale Tb v , q in a quark core given in Eq. (|B13|) . 
Thus here the contribution to damping from bulk viscos- 
ity is only a very small part of the total dissipation power, 
which is dominated by the effects from the first-order phase 
transition. 
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